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PASSIVE RANGE AND ANGLE MEASUREMENT SYSTEM AND METHOD 

FIELD OF THE INVENTION 

The present invention relates, in general, to a system and method for 
tracking a target. More particularly, the present invention relates to a method for 
passively determining range from a moving platform to an emitter, wherein the 
emitter is the target. 

BACKGROUND OF THE INVENTION 

Target acquisition and tracking are important functions in the desire 
for autonomy , in system design of a vehicle. Vehicles are often required to sense 
their environment and track targets that are crucial to their navigation profile. 
Target state estimation is required to provide or predict an accurate target state from 
a variety of sensors onboard the vehicles. 

A conventional method for producing range to a target with respect to 
a radio wave (RF) signal is to use an active sensor onboard the vehicle. Typical 
active sensors are radio or acoustic radar or laser range finder sensors. An active 
sensor transmits a signal to a target. The signal is then reflected off the target and 
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received by the active sensor. The range to the target is calculated by a processor 
onboard the vehicle based on the travel time between the active sensor transmitted 
signal and the target reflected signal. 

Passive tracking methods offer significant advantages over active 
tracking methods. Unlike radar, laser and other active tracking sensors, passive 
sensors do not emit any kind of energy. They only receive target emitted energy and 
transform the energy for measurement purposes. This characteristic makes passive 
tracking methods an ideal technique in reconnaissance and surveillance applications. 
Passive tracking methods can detect the target and, at the same time, can keep the 
detecting platform or vehicle hidden from any external detection by the target, as it 
emits no signals. 

In general, however, a passive tracking sensor cannot measure range 
or distance between the sensor and the target, as it is not based on the echoed- 
signal principle. The passive sensor, typically, offers only measurement of the target 
direction with respect to the received RF carrier in space. It is, therefore, very 
challenging to estimate accurate range information to a target from a passive sensor. 
The present invention addresses the challenge of estimating accurate range and 
angle information to a target using only passive sensors onboard a moving vehicle or 
platform. 
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SUMMARY OF THE INVENTION 

r 

To meet this and other needs, and in view of its purposes, the present 
invention provides a method of determining range from a moving platform to an 
emitter. The method includes the steps of: (a) receiving a RF signal from the 
emitter; (b) counting a number of phase reversals of the received RF signal during a 
period of time; (c) measuring a Doppler frequency during the period of time; and (d) 
determining the range to the emitter based on both the number of phase reversals 
counted in step (b) and the Doppler frequency measured in step (c). Step (b) 
includes counting the number of phase reversals of the received RF signal during the 
period of time the moving platform traverses a distance. The method may further 
include the step of: (e) measuring the distance traversed by the moving platform 
during the period of time. Step (d) may include determining the range to the emitter 
based on the number of phase reversals counted in step (b), the Doppler frequency 
measured in step (c) and the distance measured in step (e). Measuring the distance 
may include obtaining geographic position data at each end of the distance traversed 
by the moving platform, using either an inertial navigation system (INS), a Global 
Positioning System (GPS), or a combination of an INS and GPS, 

Another embodiment of the invention provides a method of 
determining range from a moving platform to an emitter. The method includes the 
steps of: (a) receiving a RF signal from the emitter during a period of time the 
moving platform traverses a distance, the distance denoted by b; (b) determining a 
carrier wavelength, X, of the RF signal; (c) counting a number of phase reversals of 
the received RF signal during the period of time, the number denoted by N; (d) 
determining a range differential, AR, between the moving platform and the emitter 
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during the period of time, in which AR=NA.; (e) measuring a Doppler frequency, fd, 
during the period of time; and (f) determining the range to the emitter based on the 
distance b, the range differential AR and the Doppler frequency fd. 

Yet another embodiment of the invention provides an apparatus, 
installed onboard a moving platform, for determining range from the moving 
platform to an emitter. The apparatus includes a receiver for receiving a RF signal 
from the emitter, an analog to digital converter (ADC) for converting the received RF 
signal into a digital signal, a memory for storing the digital signal provided by the 
ADC, a processor coupled to the memory for extracting the stored digital signal, and 

(a) counting a number of phase reversals of the digital signal during a period of time, 

(b) measuring a Doppler frequency during the period of time, and (c) determining 
the range to the emitter using both the counted number of phase reversals and the 
measured Doppler frequency. 

The apparatus may also include a GPS receiver coupled to the 
processor for obtaining geographic position of the moving platform, and the 
processor determining a distance traversed by the moving platform during the period 
of time based on the geographic position obtained from the GPS receiver. The 
apparatus may also include a mixer coupled between the receiver and the ADC for 
converting the received RF signal into an intermediate frequency (IF) signal, where 
the ADC converts the IF signal into the digital signal. The mixer may be coupled to a 
numerically controlled oscillator (NCO) for providing a coherent signal to the mixer, 
and the mixer may combine the received RF signal and the coherent signal to provide 
the IF signal. 
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It is understood that the foregoing general description and the 
following detailed description are exemplary, but are not restrictive, of the invention. 



BRIEF DESCRIPTION OF THE DRAWING 



The invention is best understood from the following detailed 
description when read in connection with the accompanying drawing. Included in the 
drawing are the following figures: 

FIG. 1 is a diagram illustrating a geometric relationship between a 
stationary ground radar and a moving aircraft in accordance with an embodiment of 
the invention; 

FIG. 2 is a diagram of the small triangle shown in FIG. 1 forming three 
sides of a, b and d, constructed in accordance with an embodiment of the invention; 

FIG. 3 is a diagram of the equilateral triangle shown in FIG. 1 forming 
the two equal sides of Ro, constructed in accordance with an embodiment of the 
invention; 

FIG. 4 is a block diagram of the system or apparatus of the present 
invention, in accordance with an embodiment of the invention; 

FIG. 5 is a graphical representation of the received radar signal as a 
function of time, in accordance with an embodiment of the invention; 
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FIG. 6 is a graphical representation of the autocorrelation amplitude 
for the received radar signal, shown in FIG. 5, in accordance with an embodiment of 
the invention; 

FIG. 7 is a graphical representation of the Doppler frequency variation 
as a function of time for the received radar signal, shown in FIG. 5, in accordance 
with an embodiment of the invention; and 

FIG. 8 is a graphical representation of an expanded main lobe of the 
Doppler frequency signal autocorrelation function for the received radar signal, 
shown in FIG. 5, in accordance with an embodiment of the invention. 

DETAILED DESCRIPTION OF THE INVENTION 

One embodiment of the invention includes measuring the range and 
angular position of a ground based RF signal source, such as a tracking radar, with 
equipment housed on an aircraft platform. With the availability of high precision 
ranging signals from Global Positioning System (GPS) and a GPS receiver onboard 
the aircraft, the invention provides a signal processing technique to measure the 
angle and range of a pulsed-doppler/continuous wave (CW) radar. The technique 
advantageously increases the aircraft's probability to complete its assigned mission 
and enhances its survivability against the ground based radar, when that radar 
constitutes a threat. The radar position or range is determined passively by 
measuring the radar signal parameters with all equipment located aboard the 
aircraft. 
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As will be explained, the method includes a step performed by a 
surveillance receiver, namely, intercepting signals from a target tracking radar (TTR) 
and storing the signals in a digital memory. The stored signals are processed to 
extract sample to sample variations of various parameter values. These measured 
(estimated) changes of parameter values are inserted into a set of equations based 
on geometric properties of two joined triangles with a common base. The geometric 
properties of the triangle, including sides and angles, are derived from the flight 
history of the aircraft and the variations of the received signal parameter values. 

The method determines the length of one side of the triangle, formed 
by the aircraft's flight path, by using a GPS receiver that is, generally, installed in the 
aircraft to compute the geolocation of the aircraft. In the method, the GPS receiver 
measures the distance flown by the aircraft between two geologically located points. 
During the flight time, for example, from point A to point B, the method, as will be 
explained, records the number of RF signal phase reversals and the Doppler 
frequency. The recorded number of the phase reversal determines the change in the 
radar range (AR) during the time taken by the aircraft to fly the distance (b) between 
end points A and B. 

As will also be explained, the invention uses the signal Doppler 
frequency. To measure the Doppler frequency of a radar emitted signal, the 
invention includes a RF signal receiver, a digital signal storage (e.g. digital RF 
memory (DRFM)) and a digital processor. An ultra stable local oscillator (LO) such as 
a numerically controlled oscillator (NCO)) is utilized for the down frequency 
translation process. 
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It will be appreciated that there is a high stability requirement on the 
LO, because the observation time is very long. The observation time may be as short 
as 1 second, or as long as 20 seconds. The longer observation time provides a 
longer range differential and a greater accuracy of the measured Doppler frequency. 

An embodiment of the method of the invention will now be described 
by reference to FIGS. 1-3. As shown in FIG. 1, moving platform 10, which may be 
an aircraft, is moving at velocity v, at an unknown range from stationary radar 20. 
The radar is emitting a RF signal. In the example shown in FIG. 1, aircraft 10 is 
flying a straight line course at a constant velocity and altitude. Locations of aircraft 
10 and radar 20 are defined in a tilted x-y plane with a reference x-axis attached to 
the aircraft. The radar is located at y=y 0 distance from the x-axis, coincident with 
the line flown by the aircraft. The range between the radar and the aircraft, as a 
function of time, may be defined as: 

r(t)=[y 0 2 +(v*t) 2 ] 1/2 

where v is the velocity of the aircraft, and t is a time variable. 

The method of the invention to locate the radar may be initiated at an 
arbitrary time. For convenience, the time variable t is set to 0. At t=0, the radar 
signal received at the receiver antenna (FIG. 4) is translated down in frequency to a 
convenient intermediate frequency (IF). The aircraft velocity includes a component 
parallel to the radar antenna beam axis (shown as a component in a direction R). 
The range R changes, and causes RF phase shifts in the received signal. By tracking 
the RF phase and counting the number of times the phase goes through 0 in T 
seconds, the change in range (shown as a) is obtained based on the number of phase 
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states N2n. As the radar RF carrier is known (for example, through electronic 
intelligence (ELINT)), the range differential becomes AR=NX. Assuming that the 
range at t=0 is Ro, then at t=T, the range increases to Ro+AR. 

The manner in which Ro is calculated will now be described. In the 
time duration T, the aircraft moves a distance b=v*T meters. At the end of the 
observation period T, the length of two (2) line segments, shown as vectors a and b 
in FIG. 1, may be determined. The b vector has the same direction as v. A direction 
is needed for the a segment, or AR segment. The angle a where AR intersects b, may 
be found based on the received signal Doppler frequency. After determining a, the 
third side (d) may be calculated based on the law of cosines formula (side, angle, 
side). The d side may then be used as a baseline for construction of an equilateral 
triangle with two long legs Ro. The equilateral triangle, not to scale, is shown in FIG. 
1. The equilateral triangle is generally designated as 35, and the small triangle is 
generally designated as 30. 

Referring next to FIGS. 2 and 3, there is shown small triangle 30 in 
FIG. 2 and equilateral triangle 35 in FIG. 3. The construction of small triangle 30 
begins with the measurements of b and AR. To find the direction of AR an angle is 
needed. The angle formed by lines b and AR (a=AR) is found from measurements of 
the Doppler frequency of the intercepted radar signal. The Doppler frequency is 
defined as: 

f d = 2v*fc/c 
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where v is the aircraft velocity, f c is the radar RF frequency carrier, and c is the speed 
of light. Since the Doppler frequency of the radar signal, measured by the receiver 
aboard the aircraft, is imprinted with one way changes in range, the Doppler 
frequency may be redefined as follows: 

f d = range-rate*fc/c 

The range rate is obtained by differentiation of the range r(t) 
expression. The analytical expression for the range rate is as follows: 

dr(t)/dt=v 2 *t/r(t) 

The Doppler frequency measured aboard the aircraft is not constant, 
and varies with time. Incorporating the geometric relationship into the expression of 
the Doppler frequency, provides the following: 

f d =v*cos(a)*fc/c 

This validates the method of the invention for determining the angular 
position of the radar from measurements of the signal Doppler frequency. Solving 
for a the following expression is obtained: 

a=Acos(f d *c/f c *v) 

where Acos is arc-cosine. 

In this manner, the angular position of the source of the RF signal 
(radar) is found. With the angle a known, the length of the third leg of the small 
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triangle bounded by b and a may also be calculated. By applying the law of cosines 
formula, the third side may be found, as follows: 

d = (a 2 +b 2 -2a*b*cos(a)) 0 ' 5 

The completed triangle Is shown in FIG. 3. With the known lengths of 
the three sides of the triangle, the angle (p) between the d and a sides may be 
calculated as follows: 

p=Acos[(a 2 +d 2 -b 2 )/2*d*a] 

where Acos is arc-cosine. 

The long legs of equilateral triangle 35 may be found next. As shown 
in FIG. 3, the angle between the base d and the long side Ro is v|/=7i-p. To simplify 
the computation of the range R, base d of the equilateral triangle may be divided into 
two halfs and a dividing line (shown as a dashed line in FIG. 1) may be drawn from 
the center of base d to the apex of the triangle, at the location of the radar. The 
dividing line forms two right triangles, with Ro as the hypothenuse. 

Finally, the method computes the range, as follows: 

Ro=d/2sin(s), and 
R=d/2*sin(e)+AR 

where e is the angle of the right triangle at the radar end. 
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A numerical example using the above described method will now be 
provided. Assume that an aircraft is flying a straight and level course at a velocity of 
200 meters per second and a radar is located at a cross range of 50km (y 0 =50km). 
The radar employs coherent signal processing techniques, operates at a frequency 
carrier of 10GHz, and has sufficient RF signal power, so that the aircraft receiver may 
detect it. An observation time (T) is set to 20 seconds. 

Twenty second observation time implies that side b is 4000 meters 
long. During the observation time of 20 seconds, 35933(2*0 phase states are 
recorded. At 10GHz frequency, the wavelength is .03 m., and the range differential 
(AR) is 1,078m. (a=l,078m.). The range differential is a vector with magnitude AR 
and has a direction to be determined. 

The direction of the range differential is derived from the one-way 
Doppler frequency. As produced by a computer model, the Doppler frequency is 
2032 Hz. Dividing the 2032 Doppler frequency by the product of the radar signal 
carrier frequency and the aircraft velocity (v*f c ), and then multiplying the result by 
the speed of light (c), the angle a is found to be 72.255degrees. It is calculated by 

a=Acos(2032*3*10 8 /200*10 10 ). 

The radar's angular position (line of sight, LOS) is 72.255 degrees with 
respect to a line parallel with the aircraft velocity vector. With the calculation of this 
angle, the direction of the range differential (length) is also found. With the two 
vectors known, the magnitude of the third side (d) is calculated using the law of 
cosines formula, as follows: 
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d = [16+l. 162-2*4*1. 078*cos(72.254)]°' 5 =3. 812km. 

To construct the equilateral triangle on base d, the angle between a 
and d is computed, using the law of cosines formula, as follows: 

(3=Acos[(14. 531 + 1. 162-16)/2*3. 812*1. 079]=92.1245degrees. 

The two equal angles of the equilateral triangle are y^-p, v|/=87.8755 
degrees. From the mid point of base d, a line normal to the base is drawn, forming 
two right triangles, as shown in FIG. 3. The second angle of the triangle is e=7i/2-v|/, 
where e=2.1245 degrees. The length of the hypothenus is Ro =d/2*sin(s)=51. 282km. 
Finally, the range R=Ro+AR may be found, resulting in R=51. 415+1. 078=52. 493km. 

The invention may be used, for example, in military type airborne 
platforms, such as tactical aircraft (fighter, reconnaissance), UAV, helicopters, etc. It 
will be appreciated that the method of the invention may be implemented using 
conventional design receivers, for example, advanced technology early warning 
(ATEW) receivers that are installed aboard many military airborne platforms. 

It will be understood that the invention may measure range and angle 
of pulsed Doppler/continuous wave (PD/CW) radars only. This class of target 
tracking radars employs coherent signals, which the ATEW receiver, for example, 
may accurately measure with the signal parameter values needed. Two parameter 
values that the invention requires are range differential and signal Doppler 
frequency. Unless the received signal is phase coherent, these two parameters 
cannot be measured. 
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Referring next to FIG. 4, there is shown an embodiment of a system 
for determining range to radar 20, the system generally designated as 40. As 
shown, system 40 includes GPS receiver 42 coupled between GPS antenna 41 and 
processor 43. Also included are electronic counter-measures (ECM) receiver 45 
coupled between ECM antenna 46 and frequency translator, or mixer 47. Mixer 47 
combines the radar RF signal from ECM receiver 45 with an oscillation signal from 
NCO 48 to produce an intermediate frequency (IF) signal on line 47a. 

The IF signal on line 47a is provided to analog-to-digital converter 
(ADC) 49 to produce a digital IF signal, which is stored in digital RF memory (DRFM) 
50 and controlled by process controller 44. Processor 51 receives the stored digital 
IF signal and computes N and the Doppler frequency. The Doppler frequency is 
computed using an autocorrelator, as described below. Processor 52 computes the 
three sides of the small triangle (b, a and d) and the length (Ro) of the two sides of 
the equilateral triangle to determine the range R to the radar. 

It will be appreciated that system 40 is installed onboard aircraft 10. 
Although shown in FIG. 4 as three separate processors, it will also be appreciated 
that processors 43, 51 and 52 may be a single processor. 

As an option, processor 43 may receive navigation data from inertial 
measuring unit (IMU)/inertial navigation system (INS) 53. The navigation data 
provided by IMU/INS 53 may be optimally integrated with the navigation data 
provided by GPS receiver 42. 

Still referring to FIG. 4, during the observation time, T, system 40 may 
measure the number of RF cycles of radar 20, the ground distance flown by the 
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aircraft (VT), and the signal Doppler frequency (fd). In the embodiment shown in 
FIG. 4, the radar signals intercepted by onboard receiver 45 are first down converted 
to an IF signal, digitized and stored in memory 50. After accumulation of signals 
over time T, the data are read out from the digital storage and processed by 
processor 51 to extract the number of cycles (N) that the signal has changed, due to 
range contraction or expansion, and the radar signal one-way Doppler frequency. 
The extraction process performed by processor 51 is further discussed below. 

One of the difficulties, in the radar location method of the invention, is 
the variable nature of the Doppler frequency and its phase instability. Changes in 
aircraft- radar geometry, causes progressive changes in the signal parameter values. 

An important parameter that processor 51 is calculating depends on 
the microscopic changes between zero crossings of the received radar signal. To 
quantify these changes, a brief analytical review of the signal is presented below. 

A radar transmitted signal waveform is of the form Vt(t)=cos(27if c t). 
The signal waveform is received by the aircraft receiver. The range between radar 
20 and aircraft 10 is expressed as R(t) = Ro+Rr(t)*t, where Rr(t) is the range rate. 
The signal waveform received at the aircraft receiver, referenced to the phase of a 
stable local oscillator, is defined as follows: 

Vr(t)=cos{t[27tf c -Rr(t)27cf c /c]-Ro27if c /c} 

where the first term is the Doppler frequency, and the second term is phase, fJc=l/\ 
and <))=Ro27iA. 
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As the range length expands, the phase of the received signal is 
varied. The phase represents the number of equivalent wavelengths to the range Ro, 
and it is a large number. Thus, the phase change is substantially equivalent to a 
distance and may be extracted by saturating the intercepted signal amplitude and 
counting the number of zero crossings. By detecting the intercepted radar signal and 
mixing it with a locally generated reference, the radar signal, as intercepted, may be 
determined. 

Extracting the signal Doppler frequency is more complex. The Doppler 
frequency under most operational conditions varies with time. If the aircraft flies a 
straight path, the Doppler frequency during short time segments varies linearly up or 
down. At the closest range, there is reversal of direction, from positive to negative. 
As one option, the Doppler frequency variation may be determined by taking a short 
segment of the signal, close to the end of T, counting the number of zero crossings 
and determining its average. The averaged Doppler frequency may be used in the 
computation of the angle a. 

As another option, the signal autocorrelation function may be 
computed. The width of the main lobe of the autocorrelation function quantifies the 
Doppler frequency signal bandwidth. This option is useful when the observation time 
is long (in seconds) and the Doppler bandwidth is in kiloHertz. 

To extract the desired radar signal parameters, a locally generated 
signal frequency is used. The incoming signal is translated down to a practical IF, 
digitized and stored in a memory, such as a DRFM. After a predetermined size 
memory is filled, the string of data is extracted and autocorrelation is performed. 
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The width of the main lobe computed autocorrelation function is proportional to the 
inverse of the time span between the 3dB points. These signal properties are 
illustrated in FIGS. 5-8 and discussed below. 

Assuming that the DRFM output signal is time compressed (eliminating 
long empty memory locations), the reconstructed output signal resembles a chirp 
signal (linear frequency modulated time waveform). The analog form of this output 
signal is shown in FIG. 5, as generated by a computer model. For simplicity, the 
observation time begins at the closest range to the radar. The signal Doppler 
frequency starts at zero (plus IF). At 50km range and 200 mps aircraft velocity, the 
Doppler frequency varies from 0 to 266 Hz in 10 seconds time interval. The output 
signal frequency is divided by 10 and graphically shown in FIG. 5. 

The autocorrelation function of this output signal, at the divided 
frequency (fd/10), is shown in FIG. 6. The Doppler frequency as a function of time 
(10 sec.) is shown in FIG. 7. The expanded main lobe of the Doppler frequency 
signal autocorrelation function is shown in FIG. 8. As shown, the 3dB power points 
of the main lobe are separated in time by 1/Afd. 

Although illustrated and described herein with reference to certain 
specific embodiments, the present invention is nevertheless not intended to be 
limited to the details shown. Rather, various modifications may be made in the 
details within the scope and range of equivalents of the claims and without departing 
from the spirit of the invention. 



